Molecular Hydrogen Kinematics in Cepheus A 
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ABSTRACT 

We present the radial velocity structure of the molecular hydrogen outflows associated to 
the star forming region Cepheus A. This structure is derived from the dopplcr shift of the H2 
v=l-0 S{1) emission line obtained by Fabry-Perot spectroscopy. The East and West regions of 
emission, called Cep A (E) and Cep A (W), show radial velocities in the range of -20 to km s^^ 
with respect to the molecular cloud. Cep A (W) shows an increasing velocity with position offset 
from the core indicating the existence of a possible accelerating mechanism. Cep A (E) has 
an almost constant mean radial velocity of -18 km s~^ along the region although with a large 
dispersion in velocity, indicating the possibility of a turbulent outflow. A detailed analysis of 
the Cep A (E) region shows evidence for the presence of a Mach disk on that outflow. Also, we 
argue that the presence of a velocity gradient in Cep A (W) is indicative of a C-shock in this 
region. Following Riera et al. (2003), we analyzed the data using wavelet analysis to study the 
line width and central radial velocity distributions. We found that both outflows have complex 
spatial and velocity structure characteristic of a turbulent flow. 



Subject headings: ISM: individual (Cepheus A) - ISM: jets and outflows - 
- ISM: molecules - ISM: Turbulence - infrared: ISM 



ISM: kinematics and dynamics 



1. Introduction 

Cepheus A is the densest core within the 
Cepheus 0B3 molecular cloud complex (Sargent 
1977) and a massive star forming region. It con- 
tains a deeply embedded infrared source which 
generates a total luminosity of ~ 2.4 x IO'^Lq 
(Koppenaal et al. 1979). 



^ Based on observations made at the 2.1 m telescope of 
tlie Observatorio Astronomico Nacional at San Pedro 
Martir, B.C., Mexico. 

^Observatorio Astronomico Nacional, San Pedro Martir, 
B.C., Mexico. 



Two main regions of ionized and molecular gas 
about 2 arc minutes apart and oriented roughly 
in the east-west direction, have been detected, 
Cepheus A East (Bally & Lane 1982) and Cepheus 
A West (Simon & Joyce 1983; Garay et al. 1996). 
The first molecular hydrogen map of both regions 
using Fabry-Perot spectroscopy was presented by 
Doyon & Nadeau (1988). Later, Hartigan et 
al. (1986) obtained high resolution images in the 
v=l~0 S{1) line emission of H2 from Cepheus A 
West. The two Cepheus A regions with molecular 
hydrogen emission show quite different composi- 
tions. 
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The eastern region (hereon Ccp A (E)) hosts 
one of the first detected CO bipolar molecular out- 
flows (Rodriguez, Ho, & Moran 1980). High reso- 
lution observations show a more complex outflow 
of a quadrupole nature (Torrelles et al. 1993). 
Torrelles et al. (1993) suggested that the source 
Ccp A East:HW2 (Hughes & Woutcrloot 1984) 
is powering the PA=45° outflow but it is not 
clear if the powering source of the PA=115° out- 
flow is Ccp A East:HW3 or another source (L. 
F. Rodriguez, personal communication). Obser- 
vations of ^^CO, CS, and CSO give evidence of 
multiple episodes of outflow activity (Narayanan 
& Walker 1996). Observations by Codella et 
al. (2003) of H2S and SO2 confirm the pres- 
ence of multiple outflows. Highly variable H2O 
and OH masers, commonly associated to young 
stellar objects, are surrounded by very dense NH3 
condensations that probably redirect the outflow 
into a quadrupole structure (Torrelles et al. 1993; 
Narayanan & Walker 1996). 

The western region (hereon Cep A (W)) con- 
tains several radio continuum sources at 3 cm 
(Garay et al. 1996) and Hartigan et al. (1986) 
identified a region of several Herbig-Haro objects 
known as HH 168 (GGD 37) with large radial 
velocities and line widths. A bipolar outflow of 
CO with overlapping red- and blue-shifted lobes 
is associated to this region (Bally & Lane 1982; 
Narayanan & Walker 1996). It should be pointed 
out that the energy source of Cep A (W) remains 
elusive (Raines et al. 2000; Garay et al. 1996; 
Torrelles et al. 1993; Hartigan & Lada 1985). 

Although the two regions Cep A (E) and 
Cep A (W) could constitute a single large struc- 
ture outflow, several authors have presented ev- 
idence which suggests that Cep A (W) may be 
an independent region of activity, distinct from 
Cep A (E) (Raines et al. 2000; Garay et al. 
1996; Hartigan & Lada 1985). 

In this paper we present the radial velocity 
structure of the Cepheus A molecular hydro- 
gen outflows obtained from the H2 v=l-0 S{1) 
doppler shifted emission line at 2.122 /im mea- 
sured by scanning Fabry-Perot spectroscopy. Due 
to the complexity found in the velocity structures, 
we decided to study the kinematics by using an 
asymmetric wavelet analysis following Riera et 
al. (2003), who used this method to study Ha 
Fabry-Perot observations of the HH 100 jet. 



Our results show that the two regions represent 
turbulent H2 outflows with significant differences 
from a kinematic point of view. A detailed analy- 
sis of the Cep A (E) region provides evidence for 
the presence of a Mach disk near the tip of the 
outflow. 

In § 2 we described the observations. From this 
data, in § 3 we generate a doppler shift H2 image, 
radial velocity, velocity gradient and line width 
maps, and study the flux-velocity diagrams (Salas 
& Cruz-Gonzalez 2002). By using the asymmetric 
wavelet transform, the clumpy structures of both 
regions of Cepheus A are kinetically analyzed and 
discussed in § 4. The conclusions are then sum- 
marized and presented in § 5. 

2. Observations 

In October 5, 1998, we observed the Cepheus A 
region with the 2.1 m telescope of the Observatorio 
Astronomico Nacional at San Pedro Martir, B.C. 
in Mexico. 

The measurements were obtained with the 
CAMILA near-infrared camera/spectrograph (Cruz- 
Gonzalez ct al. 1994) with the addition of a cooled 
tunable Fabry-Perot interferometer (located in the 
collimated beam of the cooled optical bench), and 
a 2.12 /im interference filter. A detailed descrip- 
tion for the infrared scanning Fabry-Perot instru- 
mental setup is presented by Salas et al. (1999). 

The Fabry-Perot has a spectral resolution of 24 
km s~^ and to restrict the spectral range for the 
v=l-0 5(1) H2 line emission, an interference fil- 
ter (2.122 /xm with AA= 0.02 /im ) was used. 
The bandwidth of this filter allows 11 orders of 
interference. Only one of these orders contain the 
2.122 fim line and the remaining orders contribute 
to the observed continuum. The spatial resolution 
of the instrumental array is ~ 0.86" pixel"^. The 
field of view allows to cover a 3.67' x 3.67' re- 
gion, which corresponds to 0.7 x 0.7 pc^ at the 
adopted distance of 725 pc (Johnson 1957). With 
this field of view, one set of images was required 
for the eastern portion and another for the western 
region of Cepheus A. 

Images of each region of interest were obtained 
at 26 etalon positions, corresponding to incre- 
ments of 9.82 km s~^. The observing sequence 
consists of tunning the etalon to a new position 
and imaging the source followed by a sky expo- 
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sure at an offset of 5' south from the source. 
The integration time of 60 s per frame was short 
enough to cancel the atmospheric Unes variations 
at each ctalon position, but long enough to obtain 
a good signal-to-noise ratio. Images were taken 
under photometric conditions with a FWHM of 
1.6". 

Spectral calibration was obtained by observing 
the line at 2.1332885 /xm of the Argon lamp at 
each position of the ctalon, giving a velocity un- 
certainty of 1 km s~^ in the wavelength fit. A set 
of high- and low-illumination sky flats were ob- 
tained for flat- fielding purposes. 

We reduced the data to obtain the velocity 
channel images using the software and the data re- 
duction technique described in Salas et al. (1999). 

3. Results 

3.1. H2 Velocity Maps 

Velocity channel images were individually ob- 
tained for Cep A (E) and Cep A (W) from the 
position- velocity cube data. For each pixel on the 
image we subtracted a continuum intensity level 
calculated from the median of the channels with 
no H2 emission. H2 v=l-0 S{1) line emission was 
detected in velocity channels -40 to km s~^ in 
Cep A (E) and in channels -40 to 10 km in 
Cep A (W) . The two sets of maps were pasted to- 
gether to create velocity channel maps of the com- 
plete region. Figure 1 shows five of these maps (8 
to 12) covering Local Standard of Rest (LSR) ve- 
locities from -42 km s~^ to -3 km s~^. Cep A (E) 
shows H2 emission in six separated clumps of emis- 
sion (labeled B G in Fig. 1), while Cep A (W) the 
emission may be distinguished in six regions (la- 
beled H-M). 

We crciatcd a color coded velocity image from 
the three channel velocity maps with more copious 
emission (-32.3, -22.5, and -12.7 km s~^) as blue, 
green, and red respectively. Figure 2 presents 
this color composite map. It should be noted 
that these velocities are somewhat bluer than the 
-11.2 km systemic velocity found from mil- 
limeter wavelength line observations of CO (e.g. 
Narayanan & Walker (1996)), as had already been 
noted by Doyon & Nadeau (1988). 

The velocity structure in both regions shows a 
very complex pattern. However, a slightly system- 



atic change from blue to red starting at the center 
of the image can be appreciated on the western 
outflow, that is not present in the eastern region. 
On the other hand, the Ceph A (E) region shows 
a large amount of small clumps with different ve- 
locities lying side by side. 

We have calculated the centroid radial velocity 
for each pixel by taking only 10 velocity channels 
around the peak intensity. Figure 3 presents his- 
tograms of centroid radial velocity for the two H2 
emission regions of the image in Fig. 2. Regardless 
of the difference in morphology for the two regions, 
they have similar fractional distribution of pixels 
for a given radial velocity. However, Cep A (W) 
has a wing that extends into positive radial veloc- 
ities and a small peak in the negative velocities. 

Figure 4 presents the centroid radial veloc- 
ity as a function of displacement along the right 
ascension axis for both regions Cep A (E) and 
Cep A (W). All the channels with detected emis- 
sion of H2 from Fig. 2 are shown. The radial ve- 
locity of Cep A (E) has an almost constant mean 
value of '--^-18 km along the region with a high 
dispersion around this value. Meanwhile, radial 
velocity of Cep A (W) increases its value with off- 
set position to the West. We have fitted a line 
to each position-velocity data for each region and 
found rms residual values of 8.3 and 10.2 km s~^ 
for Cep A(W) and Cep A(E), respectively. The 
smaller velocity dispersion and the large number 
of small clumps with different velocities (see Fig 2) 
in the East region indicate a more turbulent out- 
flow in Cep A(E) compared to Cep A(W). 

3.2. A Mach disk in CephA(E) 

The Cep A (E) outflow culminates in an arc 
shaped structure (labeled G in Fig. 1) that re- 
sembles a bow shock. This region is amplified in 
Fig. 5a. A bright spot can be seen in the center of 
the bow. The centroid velocities corresponding to 
this region are shown in Fig. 5c. The highest blue- 
shifted velocity of the region (-40 km s~^) corre- 
sponds to a slightly elongated region (in the direc- 
tion perpendicular to the outflow) that includes 
the bright spot and decreases toward the bow, as 
can also be seen in a position- velocity diagram in 
Fig. 5d. This kinematic behavior is expected if 
the bright spot corresponds to the Mach disk of 
the jet, where the jet material interacts with pre- 
viously swept material accumulating in front of the 
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Fig. 1. — Velocity channel maps of the H2 v=l 5(1) line emission in Ccpheus A, showing the radial 
velocity structure of the two regions Cep A (E) and Cep A (W). The velocity varies from -42.1 km s~^ in 
channel 8 (top) to -2.9 km s~^ in channel 12 {bottom), the channel width is 9.8 km s~^. The scale of the 
region is shown at the top panel. The contour intervals indicate the emission intensity in the range of 0.54 
to 1.5 counts/sec and have the same increment values of 0.2 counts/sec. The cross indicates the position of 
the 6 cm peak of HW2 at R.A. = 22h45ml7.9s and DEC. = -^62° 01' 49" (2000), see Hughes & Wouterloot 
(1984). Offsets are referred to this reference position. 
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Fig. 2. — Color composite of the doppler shifted H2 hne emission t;=l-0 S{1) in Cepheus A, showing the 
two main regions Cep A (E) and Cep A (W). Colors correspond to velocities in the LSR: Blue -32.3 km s~^, 
green -22.5 km s~^, and red -12.7 km s~^. The region shown is about 5'x3', North is at the top and East is 
to the left. 
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Fig. 3. — Fraction of pixels with a given value of the centroid radial velocity (LSR) for the Cep A (E) {dashed 
line) and Cep A (W) {solid line) regions. The size of the bins is 2 km s"^. 



jet and behind the bow shock. 

A few cases arc known where a Mach disk is 
observed in H2. The case of Kumar et al. (2002) 
at the Nl outflow in S233IR (Porras et al. 2000), 
where a flattened structure is seen in high defini- 
tion H2 images, although it is not supported spec- 
troscopically. In HH 7 Khanzadyan et al. (2003) 
detect a flattened [Fell] structure that coincides 
with a blue-shifted knot in H2. This case is simi- 
lar to the present case for Cep A (E). Blue-shifted 
velocity is expected as the jet material flows in all 
directions upon interacting at the Mach disk, in 
particular toward the present direction of the ob- 
server, nearly perpendicular to the outflow axis. 
As noted by Kumar et al. (2002) a detection of 
a Mach disk has rather interesting implications: 
1) the jet would be partially molecular; 2) the ve- 
locity of the jet must be small enough to prevent 
dissociation of H2 in the Mach disk; and 3) the jet 
should be heavy. In this case wc can estimate the 
velocity of the molecular jet from the maximum 
observed velocity (-40 km s~^) minus the rest ve- 
locity of the molecular cloud (Vb=-11.3 km s~^) to 
be around 29 km s~^. The distance from the Mach 
disk to the bow-shock apex is 5" to 8" (~0.016 to 
0.025 pc at D=725 pc). 



3.3. Cep A (W) velocity gradients 

The western H2 outflow in Cep A displays a 

series of wide arcs reminiscent of thin sections of 
shells. The outflow has been described as a hot 
bubble (Hartigan, Morse, & Bally 2000) that 
drives C-shocks into the surrounding medium. 
In some shock fronts along this outflow, they 
observed that H2 emission leads the optical 
[SIIJA6717 which in turn leads Ha. This is taken 
as evidence in favor of a C-shock that slowly ac- 
celerates and heats the ambient medium ahead 
of the shock. We find further evidence of this 
from the H2 kinematics. A velocity gradient from 
higher to lower velocities is observed in some of 
the individual arc structures, as is the case for 
the one labeled I in Fig. 1. The spatial map of 
the centroid velocities in Cep A (W) is shown in 
Fig. 6. A smooth velocity gradient, covering from 
-36 km s~^ to -8 km s^^, is observed in the SW 
border of the arc, in the direction indicated by 
the arrow. That is, a velocity gradient going from 
blue-shifted to closer to the rest velocity of the 
molecular cloud (-11.2 km s~^) in a region of 17", 
or 0.06 pc at D=725 pc. We regard this as the 
kinematic evidence of a C-shock, as the slow ac- 
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Fig. 4. — Diagrams of detected H2 v=l-0 S{1) line emission in all velocity channels(top) and the velocity- 
position diagrams {hottom) for Ccp A (E) (left) and Cap A (W) {right). Offset values are referred to the 
position of HW2 (see Fig. 1) considered as position (0,0). 
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Fig. 5. — Mach disk in Cep A (E). (a) Amplified H2 image of region G showing the edge of the outflow 
presented in panel (b). (c) shows the centroid velocity of the bow-shock region G in gray scale, going from 
—40 km s~^ (black) to km in white. Overlaid contours show the intensity map. (d) shows a position- 
velocity diagram of a cut shown in (c) in the direction indicated by the arrow. The solid line corresponds 
to velocity —{V — Vo), Vq = —11.3 km s^^, and the dotted hne is intensity. 
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celeration of material ahead of the shock is quite 
evident. 

3.4. Flux— velocity Relation 

Wc have calculated the flux-velocity diagrams 
separately for the Cep A (E) and Cep A (W) 
outflows, as described in Salas & Cruz-Gonzalez 
(2002). For every pixel with signal above the de- 
tection threshold, we added the fluxes of all the 
pixels with centroid velocities in bins of observed 
centroid velocity l?;,,^,, — Urest |, with respect to the 
rest velocity of the region (-11.2 km s~^). The 
flux-velocity diagrams so obtained are shown in 
Fig. 7 and 8. 

As shown by Salas & Cruz-Gonzalez, this pro- 
cedure gives similar flux- velocity relations for a va- 
riety of outflows, consisting in a flat spectrum for 
low velocities, followed by a power-law decrease 
above a certain break-velocity. The power-law in- 
dex is very similar for different outflows, as is the 
case for Cep A (E) and Cep A (W) where it is 
-2.6±0.3 and -2.7±0.9 respectively (solid lines in 
figures 7 and 8). The break velocity however, is a 
little different. The logarithm of Vf,reafc (inkms~^) 
takes values of 0.95±0.07 and 1.13±0.13 respec- 
tively, a difference of around 2a which suggests 
that Vbreak may be larger for Cep A (W). As was 
discussed in Salas & Cruz-Gonzalez (2002), out- 
flows of different lengths (Z) show break-velocities 
varying as Ubreafc »: a result that is taken to 
imply an evolutionary effect, similar to the case 
of CO outflows (Yu et al. 1999). However, in 
the case of Cep A (E) and Cep A (W) the out- 
flow length is very similar, as might be the outflow 
age. Salas & Cruz-Gonzalez also argue that other 
causes for a difference in break velocities could be 
the amount of turbulence in the outflow, which 
might be the case for Cep A as is mentioned at 
the end of §3.1. We will next explore this possi- 
bility through the use of a wavelet analysis. 

4. Wavelet Analysis of the H2 emission 
from Cep A 

4.1. Description 

Both H2 emission regions in Cep A show a 
rather complex velocity-position structure (see 
Fig. 2). We have carried out a wavelet analysis 
in an attempt to understand the relation between 



sizes of clumps, velocity, and velocity dispersion 
as a function of position along the outflows. 

The wavelet transform analysis has been used 
before by Riera et al. (2003) to analyze Ha ob- 
servations using a Fabry-Perot of the complex out- 
flow in the HH 100 jet. Although, their procedure 
is equivalent to the one followed by Gill & Hen- 
riksen (1996) to study turbulence in molecular 
clouds, we have followed the analysis by Riera et 
al. (2003) because it is more relevant to the study 
of outflows. 

The two-dimensional wavelet transform analy- 
sis is used to obtain the positions of all H2 emission 
clumps as well as their characteristic sizes. The 
sizes of these clumps are determined from total H2 
flux images created by adding the continuum sub- 
tracted intensities of all the channel maps. Once 
obtained the positions and characteristic sizes of 
the H2 structures of Cep A, the spatial averages 
(over the characteristic sizes of the structures) of 
the line center velocity and of the line widths are 
calculated to study the turbulence in the two re- 
gions of Cep A. 

Our case is similar to Riera et al. (2003) in 
the sense that Cep A regions have an axial sym- 
metry and that the kinematics properties seem to 
depend on the offset along this axis. Then, it is 
reasonable to expect that the characteristic size 
may be different along and across the regions. To 
exploit this symmetry, we rotated the H2 image 
so the long part of both regions, the main axis, 
is parallel to the x coordinate defined along the 
east-west direction. On this rotated images, we 
carried out an anisotropic wavelet decomposition 
using two bases of different sizes, one along and 
one across the main axis. 

We adopted a basis of "Mexican hat" wavelets 
of the form 

gir;a,,ay) = Ci2-r^)e-^'/\ (1) 

where r = [(x/ax)^ + (y /ay)"^]^^^; ax and ay are 
the scale lengths of the wavelets along the x- and 
y-axes, respectively, and C = {a^ + ay)^^^^. This 
is a very common used wavelet, but we also choose 
to use it because simplifies the detection of inten- 
sity peaks, approaches better to the shape of the 
intensity peaks, and behaves well under FFT cal- 
culations. 

To compute the wavelet transform we have to 
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Fig. 6. — Velocity gradients in Cep A (W). The gray levels show velocity in the range -40 km s ^ (black) to 
-4 km s~^ (white). 
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Fig. 7. — Flux-velocity relation of H2 t;=l-0 S{1) line emission in Cep A (E). The long- dashed line shows 
the flat low velocity behavior up to a break velocity where the solid line indicates a power law decrease. The 
dashed lines show the range of possible values of the index (-2.6±0.3) and of log{vbreak) (0.95±0.07). 
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Fig. 8. Same as in Fig. 7, but for the region Cep A (W). The range of possible values of the power law 
index is -2.7±0.9 and of log{vbreak) is 1.13±0.13. 
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calculate the convolutions 

Ta^,ay{x,y) = j j I{x\y')g{r' ■,ax,ay)dx' dy' , 

(2) 

for each pair {ax, ay), where r' = {[{x' — x)/ax]'^ + 
[{y' ~ y)/0'y]^Y^^j I{x,y) is the intensity at pixel 
position {x,y), and g{r';ax,ay) is give by eq. (1). 
These convolutions arc calculated by using a Fast 
Fourier Transform (FFT) algorithm (Press et al. 
1992). 

The wavelet transformed images Ta^,ay{x,y) 
correspond to smoothed versions of the intensity 
of the H2 image. We use this images to find the 
sizes of the structiucs in the H2 regions of Ccph A. 
First, on the transformed image with = = 1, 
we fixed the position of x and find all the values 
of y where Ta^,ay{x, y) has a loc;al maximum. Sev- 
eral maxima may be found for each position x that 
correspond to different structures observed across 
the regions. The maxima found with ax = dy = i 
will correspond also to the local maxima of I{x, y). 

For each pair where I{x, y) has a max- 

imum we determine {ax,ay), in the and ay 
space, where the wavelet transform has a local 
maximum, a^ and ay will be the characteristic size 
of the c;luinp with a maximum intensity at {x,y). 
The {ax,ay) space is search in such way that we 
first identify the size of the smaller clumps and 
then the bigger ones. This progressive selection 
allow us avoid to choose clumps that overlap with 
its neighbors. Naturally, the biggest clumps will 
have an structure similar to the whole region. 

4.2. Size of the H2 Clumps 

The results obtained with the process described 
above are shown in Fig. 9. This figure shows 
the two images of molecular hydrogen emission 
of Cep A, which have been rotated by 19° so the 
longways dimension of the regions are more or less 
parallel to the x-axis. In the case of Cep A (E), 
the X axis has been also inverted (West is to the 
right) so the x values in both panels in Fig. 9, al- 
though arbitrary in origin, are an estimate of the 
offset from the central region between the two re- 
gions. It has to be clarified that the x coordinates 
are values from independent images so there is re- 
ally no correlation between them. However, the 
span in the vertical and horizontal axis are kept 
the same in both panels of Fig. 9 to ease the size 



comparison for each region. Six large structures 
(B to G) were identified in Cep A (E) and six (H 
to M) for Cep A (W) (see Fig. 1). The spatial 
limits of these regions are shown in Table 1. 

The H2 intensity maps have then been con- 
volved with a set of wavelets g{r;ax,ay) with 
1 ^ Oa: < 30 pixels and 1 < Oy < 30 pix- 
els for Cep A (E) and 1 < a^; < 35 pixels and 
1 < a?/ < 35 pixels for Cep A (W), with a reso- 
lution of 0.853 arc-sec per pixel in the .r-direction 
and 0.848 arc-sec per pixel in the y-direction. The 
upper limiting values for ax and ay were selected 
to allow very few overlapping in the size of the 
regions in adjacent peaks. 

The values of the peaks yk of the wavelet trans- 
form obtained for each position x along the region 
are shown in the bottom panel of Fig. 9. Also 
shown as error bars are the values of ax,k and ay^k 
corresponding to each peak, which give an esti- 
mate of the characteristic sizes along and across 
the observed structures. 

To quantify the observed broadening of the re- 
gion, for each position x wc compute the weighted 
mean of the y-spatial scale perpendicular to the 
longest axis of the regions: 



Z^k -'■ax,k,ay,k{x,yk) 

Figure 10 shows (ay) as a function of position x 
along the two Cepheus A regions. In both outfiows 
we notice the presence of small (2 pixels) and large 
(25 pixels) structures located with no discernible 
order along the axial position. 

The values of function of ax,k are 

shown in Fig. 11 Cep A (E) and Cep A (W). These 
graphs represent a measure of the symmetry of the 
structures. Isotropic structures are expected to at- 
tain ay^k = o,y,k and should thus be located on a 
line of unitary slope. The values of the slope a for 
each region are also given in Table 1. We note that 
all the values are in general less than one, indicat- 
ing structures that are longer than wider. The 
western outflow, however, shows relatively wider 
structures than Cep A (E), a property which has 
also been qualitatively observed. It may be possi- 
ble that the individual structures follow the pat- 
tern of each region as a whole, with the eastern re- 
gion being longer than wider while the western re- 
gion is the opposite. However, the nature of these 
anisotropics it is not clear. 
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Fig. 9. — Position and characteristic sizes of the structures {bottom) and contour plots of the H2 image (top) 
obtained from the velocity channels for Cep A (E) {left) and Cep A (W) {right). Contour levels shown from 
0.0 to 1.5 counts/sec in 0.1 counts/sec increments for both plots. The plus sign shows the positions of the 
maximal yk of the wavelet transform obtained for the different values x along the region. The characteristic 
sizes ax,k and ay^k of these maximal are shown as error bars centered on the positions of the maximal. We 
have separately rotated the map of each region by 19° so they are more or less parallel to the a;-axis. The 
distances y (across the regions) and x (along the regions), in pixels units, are with respect to the rotated 
image and have an arbitrary origin. In the Cep A (E) case the x axis has been inverted (West is to the right) 
for analysis purpose (see text). 
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Fig. 10. — Values of (uy) plotted as a function of position x along the Cep A (E), top panel, and Cep A (W) 
at the bottom one. 
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Fig. 11. — Left: Values of Uy^k plotted as a function ax,k for the different regions along Cep A (E). The 
linear fits to the points are drawn as solid lines. The slope m of the linear fit is listed in Table 1. Right: 
Same but for the region Cep A (W). 
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Table 1: H2 Emission Regions in Cep A 



Region 






Vmin 


Vmax 


m 


a 


r> 
a 


TO 


QQ 

00 




iio 




-U.ulo: 


c, 




t7 -L 


129 




n 774 


140 


D 


92 


119 


133 


170 


0.702 


0.169 


E 


127 


140 


142 


158 


0.739 


0.155: 


F 


140 


153 


144 


155 


0.394 


0.644 


G 


163 


199 


132 


160 


0.768 


0.249 


H 


67 


92 


129 


144 




0.397: 


I 


92 


110 


129 


155 


0.871 


-0.089: 


J 


113 


129 


132 


147 


0.654 


0.563 


K 


130 


142 


115 


162 




0.379 


L 


142 


187 


115 


150 


0.931 


0.439 


M 


188 


206 


113 


131 


0.937 


0.376 



Col.(l).- Region B-G Cep A(E); Region H-M Cep A(W) (see Fig. 9). 

Col.(2)-(5) ■~ Xminy 3:rnax^ ymini d-ti-d Vmax the boundary coordinates (in pixels) of the region (see Fig. 9). 

Col. (6) .- m is the slope of the linear fit ay = max + b, where ax and ay are the sizes of the identified H2 structures on the 

region (see § 4.2) 

Col. (7).- a is the slope of the linear fit (Ai;2)i/2 = a{al + a\YI'^ + /3, where (Au^) V2 ig the rms velocity dispersion on the 
structure of size and ay (see § 4.4) 



4.3. Spatial distributions of the radial ve- 
locities and the line widths 

In this section we describe the spatial depen- 
dence of the kinetic properties of the two main H2 
emission regions of Cep A. 

From the cube of position- velocity data, we cal- 
culate two moments of the line profiles for each 
pixel: 



2003): 



Vr 



and 



J lydV ' 

J{r-V,:fl,dv 



(4) 



(5) 



/ lydV 

In these equations, v is the radial velocity, and ly 
the intensity at a fixed position {x, y) of succes- 
sive channel maps. The integrals are carried out 
over all the velocity channel maps. Here Vc is the 
barycenter of the line profile (i.e. the "line center" 
radial velocity), and is a second-order moment 
that reflects the width of the line profile. 

With these values of Vc and W computed for 
all positions (a;, y) on the plane of the sky, we cal- 
culate the following spatial averages (Riera et al. 



{Vo) = 



Is Vc{x',y')Iix',y')dx'dy' 
I{x',y')dx'dy' 



J,_^^W^{x',y')I{x',y')dx'dy' 
~fs nx',y')dx'dy' 



(6) 



, (7) 



[Vcix',y')~{Vc)]^Iix',y')dx'dy' 



Is I{x',y')dx'dy' 



(8) 

where I(x',y') is the II2 flux obtained from co- 
adding all the channel maps. These integrals are 
carried out over areas Sa^,ay which are ellipses 
with central positions (x, y) and major and minor 
axes, dx and tty, corresponding to all the values 
that have been identified as peaks of the wavelet 
transform. 

The value of {W'^){x, y) corresponds to the Hue 
width spatially averaged over the ellipse Sa^ay 
with a weight I{x',y'). The same weight spatial 
average is calculated for the line center velocities 
{Vc) within the ellipse, as well as the standard de- 
viation {Av^y^'^ of these velocities. 
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Figure 12 shows the hne centers, widths, and 
standard deviations as a function of position x 
along the region (all of the points at different posi- 
tions y across the region and with different and 
Uy are plotted). The line center or centroid ve- 
locity displayed in the top panels, shows different 
behaviors for Cep A (E) and Cep A (W). As had 
been previously noted, Cep A (E) has a constant 
velocity along the outflow (~-19 km s~^) while 
Cep A (W) shows a velocity gradient from about - 
21 to -2 km s~^. A large dispersion in the line cen- 
ters for Cep A (W) around positions 115 and 162 
is probably due to insufficient S/N, since a similar 
increase in the dispersion velocity (medium panel) 
is concurrent, but lacks a corresponding increment 
of the line width (bottom panel). Other features 
present in these figures seem real. Most remark- 
ably, the velocity dispersion (medium panel) for 
Cep A (E) increases monotonically with distance 
from the source, going from around 4 to 10 km s~^ 
in 150 pixels or 0.5 pc. No such increase is ob- 
served in the western outflow, where we obtain 
a constant 3±2 km s~^ velocity dispersion. The 
line widths (bottom panel) seem dominated by the 
width of the instrumental proflle of 22 km s~^. 

Hence, we conclude from this analysis that the 
eastern and western outflows seem intrinsically dif- 
ferent. While the eastern outflow shows a constant 
line center velocity and an increasing velocity dis- 
persion, the western outflow behaves otherwise, 
showing a velocity gradient and a constant veloc- 
ity dispersion with a lower value. This had been 
noticed in our qualitative analysis of the obser- 
vations. Now, if we take the velocity dispersion 
within each cell as a measure of turbulence, which 
seems reasonable, then the region Cop A (E) is 
more turbulent than Cep A (W), and it is notice- 
able that the turbulence increases with distance 
from the "central" source compared to a constant 
behavior in the western source. 

4.4. Deviations of the line center velocity 
and the size of the region 

Figure 13 shows the deviations of the line center 
velocity (velocity dispersion), averaged over sizes 
chosen from the wavelet spectrum, as a function of 
the size of each region. Such deviations of the line 
center velocity, averaged over sizes chosen from a 
wavelet spectrum, have been previously used by 
Gill & Henriksen (1996) to study turbulence in 



molecular clouds. In this figure wo can appreciate 
a large dispersion of values. However, points tend 
to clump together around certain regions of the 
diagrams. 

A closer examination of Fig. 13 for Cop A (E) 
shows a cluster of points at position (20,9) corre- 
sponding to structure labeled G (see Fig. 9), while 
the points corresponding to structure C lie just be- 
low of them. Taken separately, each one of these 
groups seem to define a line in the diagram, off- 
set vertically from one another, but of compara- 
ble slopes. The slope of the lines, in fact, is also 
similar to the slope of the lino that would be ob- 
tained by fitting a line to all the data, rn=0.26 
(indicated with a solid line in Fig. 13). The dis- 
persion is largo, but there is at least another fac- 
tor contributing to the amount of velocity disper- 
sion that has already been mentioned: the dis- 
tance from the components of the structure to the 
outfiow source. We take out this effect by sub- 
tracting this systematic linear contribution with 
the distance. Thus, in Fig. 13, points belonging 
to the C group, which are closer to the source, 
would merge together with those points of the far- 
ther G group. The dispersion would bo lower and 
the slope would have a similar value. This defines 
a relation similar to Larson's law for molecular 
clouds (Larson 1981). Larson found that molec- 
ular clouds show a suprathermal velocity disper- 
sion that correlates with the size a of the cloud as 
ay oc o°'^^. Although different values of the power 
index, a, have been mentioned in the literature, 
ranging from 0.33, for pure Kolmogorov turbu- 
lence, to 0.5, a more appropriate value in terms 
of the Virial Theorem (Goodman et al. 1998). In 
the Cop A (E) case we obtain a slope a of 0.25 
which is closer to the Kolmogorov value. 

In the western outflow it is more difficult to 
identify clumps of points, and we were not able 
to find a relation of velocity dispersion with posi- 
tion either. A least square fit of a line to the data 
gives a slope of 0.20. However, there appears to 
be two disperse clumps of points, a lower ono and 
an upper one, each defining a line with a larger 
slope a ~0.5 and parallel to each other. This lat- 
ter value would be closer to a case where Virial 
Theorem applies. 

The analysis of the individual H2 condensations 
(c.f. Fig. 1) in the Cep A (E) and Cep A (W) is 
illustrated in Fig. 14 and 15. We show the devi- 
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Fig. 12. — Average line centers {top panels), standard deviations (middle), and line widths {bottom) as a 
function of position x along the regions Ccp A (E) (left column) and Cop A (W) (right column). All points 
with different positions y across the region and with different and Uy are plotted. 
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Fig. 13. Relation of the deviations of the line center velocity, (A?;^)-^/^ (or a), averaged over regions chosen 
from the wavelet spectrum, and the size a of the region, (a^ +a^)^/-^, for the complete Cep A (E) (left) and 
Cep A {W){right) regions. The solid line is the best fit to the points with a slope a of 0.2574 for Cep A (E) 
and 0.199 for Cep A (W). 
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ations for the line center velocity, averaged over 
sizes chosen from the wavelet spectrum for each 
condensation. For each H2 condensation the val- 
ues of the best fit slope a are presented in the last 
column of Table 1. The Cep A (E) clumps yield a 
mean a of 0.21±0.21 while the Cep A (W) yield a 
slightly higher value of 0.34±0.20. However, some 
knots show poor correlations than others (see plots 
for knots B, E, H and I). Using only knots C, D, 
F and G for Cep A (E) and knots J, K, L and 
M for Cep A (W) yields an a of 0.30±0.20 and 
0.44±0.08, which support a Kolmogorov case in 
the first and a more virialized region for the west- 
ern case. 

5. Conclusions 

We have presented the velocity structure of 
the molecular hydrogen outfiows from the regions 
Cep A (E) and Cep A (W) obtained from the 
H2 v=l-0 5(1) doppler shifted fine emission at 
2.12 /im. Both the velocity channel maps and the 
integrated H2 image show a complex structure of 
12 individual clumps along two separated struc- 
tures oriented roughly in the east-west direction. 

Given the complexity of these structures, we 
have carried out an anisotropic wavelet analysis 
of the H2 image, which automatically detects the 
position and characteristic sizes (along and across 
the region axis) of the clumps. 

1. There is evidence for a Mach disk in Cep A (E). 
The efflux point is located at the center of a 
bow shock structure and we measure blue- 
shifted velocities of 22 to 28 km s~^. This 
observation implies that a molecular jet is 
driving the outflow. 

2. Cep A (W) on the other hand, is consistent 
with a hot bubble in expansion driving C- 
shocks. We presented the kinematic gradient 
of one of such shocks as an example. 

3. The H2 fiux-velocity relation is present in 
both outflows. The break velocity of the 
eastern outflow is lower than that of the 
western outflow, and we have argued that 
this is indicative of greater turbulence in 
Cep A (E). 

4. The wavelet analysis has conflrmed and 
quantified trends observed in the centroid 



velocity measurements: that the eastern 
outflow shows a constant line center velocity 
and an increasing velocity dispersion, while 
the western outflow shows a velocity gradi- 
ent and a constant velocity dispersion. The 
larger velocity dispersion and gradient in the 
eastern outflow is taken as indicative of tur- 
bulence, and allows us to conclude also that 
turbulence increases with distance. 

5. Suggestive propositions about the kind of 
turbulence present in both outflows, are ex- 
tracted from an analysis of the relation of the 
velocity dispersion as a function of the size of 
the structures (cells) identified as unities by 
the wavelet spectrum. Using only knots with 
a good correlation yields an a of 0.30±0.20 
for Cep A (E) and 0.44±0.08 for Cep A (W), 
which support a Kolmogorov case in the first 
and a more virialized region for the western 
case. 
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spectrum, (At)^)^/^, and the size of the region + fly)^^^ for regions B to G of Cep A (E). The solid line 
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